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ABSTRACT 

In Selective Harmonic Elimination-Pulse Width Modulation (SHE-PWM) method, low order harmonics are 
eliminated, while the fundamental harmonic is obtained at the desired value. In this paper. Variable Sampling Ant Colony 
Optimization (SamACO) algorithm with random initial values is proposed for solving the transcendental nonlinear 
equations known as Selective Harmonic Elimination (SHE) equations that characterize the selected harmonics in an 11- 
level inverter. The algorithm is a continuous (combinatorial) optimization algorithm that is based on the food foraging 
behavior of ants in a swarm. The dynamic exploitation and random exploration operators in the algorithm ensure both 
accuracy and convergence to global optima. Fast Fourier Transform (FFT) analysis of the synthesized voltage waveform 
reveals the complete elimination of the 5'*, 7'*, 1 1'* and 13'* harmonics as their values tend towards zero. Both 
computational and MATLAB simulation results show that the proposed method is highly efficient for elimination of the 
selected low order harmonics as well as minimization of the total harmonic distortion (THD). 

KEYWORDS: Multilevel Inverter, Samaco, Selective Harmonic Elimination (SHE), and THD 

INTRODUCTION 

A multilevel voltage source inverter is a power electronic system that synthesizes a nearly sinusoidal output 
voltage from several DC voltages using pulse width modulation (PWM) technique. Due to the smaller voltage steps in the 
output staircase waveform, multilevel inverter has many advantages over the traditional among which are: improved power 
quality, low switching losses, lower dv/dt stresses on the load, lower electro-magnetic interference (EMI) and ability to 
attain a higher voltage without the use of transformer [1]. The use of multilevel inverter is prevalent in industrial 
applications such as drives, Flexible AC Transmission Systems (FACTS), Hybrid Electric Vehicle (HEV), High Voltage 
Direct Current (HVDC) lines. 

The concept of multilevel inverters was actually developed from the idea of step approximation of sinusoid [2]. 
Basically, there are three main multilevel topologies. These are Diode-Clamped Multilevel Inverter [3], Capacitor- 
Clamped Multilevel Inverter [4], and Cascaded H-bridge Multilevel Inverter with separate DC sources [5]. Varieties of 
each topology as well as hybrid of the fundamental topologies such as Generalized P2 Converter, Mixed-Level Hybrid 
Converter, Asymmetric Hybrid converter have been developed but with the same underlying principle^, 7]. 

Several pulse width modulation techniques used in conventional two-level inverter have been modified and 
deployed in multilevel inverters. These include Sinusoidal Pulse Width Modulation (SPWM), Selective Harmonic 
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Elimination (SHE) method. Space Vector Control (SVC), and Space Vector Pulse Width Modulation (SVPWM) [6, 7]. 
Selective Harmonics Elimination (SHE) method at fundamental switching frequency however, arguably gives the best 
result because of its high spectral performance and considerably reduced switching loss. The main challenge associated 
with the SHE method is how to obtain the analytical solutions of the transcendental nonlinear equations that characterize 
the selected harmonics in multilevel inverter. 

The convergences of the classical derivative-dependent solvers such as Newton Raphson method [8] are highly 
sensitive to arbitrarily chosen initial values of the solution. Another approach uses Walsh functions [9] where solving 
linear equations, instead of non-linear transcendental equations, optimizes the switching angle. The method results in a set 
of algebraic matrix equations and the calculation of the optimal switching angles is a complex and time-consuming 
operation. Chiasson et al [10] proposed a method based on Elimination theory using resultants of polynomials to determine 
the solutions of the SHE equations. A difficulty with this approach is that as the number of levels increases, the order of the 
polynomials becomes very high, thereby making the computations of solutions of these polynomial equations very 
complex. 

Population-based Evolutionary Algorithms (EAs) such as Genetic Algorithm (GA) [11]-[13], Ant Colony System 
(ACS) [14], Bee Algorithm (BA) [15], Particle Swarm Optimization (PSO) [16] have been used to solve SHE equations. 
EAs are simple, derivative free and can be used for problems with any number of levels. They are also successful in 
locating the optimal solutions, but they are usually slow in convergence and require much computing time. 

MULTILEVEL INVERTERS 

• Cascaded H-bridge Inverter 

Among multilevel inverter topologies, cascaded H-bridge inverter requires the least number of components. Its 
modular structure as well as circuit layout flexibility makes it suitable for high voltage and high power applications. 
Cascaded H-bridge multilevel inverter is formed by connecting several single -phase H-bridge inverters in series as shown 
in Figure 1. The number of output voltage levels in a cascaded H-Bridge inverter is given by N = 2S + 1 , where S is the 
number of H-bridges per phase connected in cascade. By different combinations of the four switches 5] , S 2 , S 3 , and S 4 
shown in the Figure 1, each H-bridge switch can generate a square wave voltage waveform on the AC side. To obtain 
+ V dc , switches 5, and S 4 are turned on, whereas -V dc can be obtained by turning on switches .S', and S, . By turning on 
.S', and .S', , or 5, and .S 4 , the output voltage is zero. 




Figure 1: Configuration of an 11-Level Single-Phase Cascaded H-Bridge Multilevel Converter 
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The outputs of H-bridge switches are connected in series such that the synthesized AC voltage waveform shown 
in Figure 2 is the summation of all voltages from the cascaded H-bridge cells [8, 10]. 




Figure 2: Output Voltage Waveform of an 11-Level Inverter 
• Selective Harmonic Elimination PWM 

Generally, any periodic waveform such as the staircase waveform shown in Fig. 2 can be shown to be the 
superposition of a fundamental signal and a set of harmonic components. By applying Fourier transformation, these 
components can be extracted since the frequency of each harmonic component is an integral multiple of its fundamental 

[17]. 

Assuming a quarter wave symmetry and the equal amplitude of all DC sources, the Fourier series expansion of the 
staircase output voltage waveform shown in Fig. 2 is given by equation (1). 

Vicot) = V n (ar)sin(/7*Mf) (1) 

Where 



V„ ( a ) = cos (na k ) , for odd n (2) 

H7T 

V n (a) = 0 , for even n (3) 

In three-phase power system, the triplen harmonics in each phase need not be cancelled as they automatically 
cancel in the line-to-line voltages as a result only non-triplen odd harmonics are present in the line-to-line voltages 8] 

Combining equations (1), (2) and (3), 

4V 

v(ca) = X T n 5 — — (cos^a,)+cos^ar 2 )+... +cos(na^))sin neat) ... 



Subject to 0 < a t < a 2 < ...a s < ^ 

Where, S is the number of switching angles and n is the harmonic order. Generally, for S number of switching 
angles, one switching angle is used for the desired fundamental output voltage V and the remaining (5-1) switching angles 
are used to eliminate certain low order harmonics that dominate the Total Harmonic Distortion (THD) such that equation 
(4) becomes 

V((ot) = V, sin(<Mr) (5) 
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From equation (4), the expression for the fundamental output voltage V , in terms of the switching angles is given 
by 

4V i \ 

V, = — — [cos(a, ) + cos(a, ) + ...+ cos(a 5 ) ) 

K ' ( 6 ) 

The relation between the fundamental voltage and the maximum obtainable fundamental voltage Vi max is given by 
modulation index. The modulation index, mi, is defined as the ratio of the fundamental output voltage V , to the maximum 
obtainable fundamental voltage V lmax . The maximum fundamental voltage is obtained when all the switching angles are 
zero [8]. From equation (6), 



4 SV± 

71 



(7) 



m i 



V, 



7ft, 

4 



Hence, 



V, = m j 




for 0 < m, < 1 



( 8 ) 



To develop an 11-level cascaded multilevel inverter, five SDCSs are required. The modulation index and 
switching angles that result in the synthesis of AC waveform with the least Total Harmonic Distortion (THD) can be found 
by solving the following transcendental nonlinear equations known as SHE equations that characterize the selected 
harmonics[8], [10]: 

4V 

— — (costa, ) + cos(ar, ) + ... + cos(a 5 )) = V, 

K 

cos(5a, ) + cos(5a, ) + ... + cos(5a s ) = V 5 

cos(7ar, ) + cos(7a, ) + ...+cos(7a 5 ) = V 1 

cos(l la , ) + cos(l la, ) + ... + cos(l la 5 ) = V , , 

cos(13a, ) + cos(13a, ) + ... + cos(13a 5 ) = V I3 (9) 

In equation 9, V 5 , V 7 , V n , and V 13 are set to zero to in order to eliminate 5 th , 7 th , 11 th and 13 th harmonics 
respectively. The correct solution must satisfy the condition 

0<a, <a 2 <...<a 5 <^>2 (10) 

Equation (8) in equation (9) yields: 

cos(a, )+ cos(a, ) + ...+ cos(a 5 ) = 5m, cos(5a,) + cos(5a,) + ...+cos(5a 5 ) = 0 cos(7a,) + cos(7a,) + ...+cos(7a 5 ) = 0 
cost 13a, ) + cos(l 3a, ) + ... + cost 13a 5 ) = 0 (11) 

Generally equation (11) can be written as 
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F(a) = B(m i ) 



(12) 



The Total Harmonic Distortion (THD) is computed as shown in equation (13): 




(13) 



ANT COLONY OPTIMIZATION (ACO) 



The ant colony optimization is swarm intelligence based meta-heuristic algorithm that was inspired by the food 



the non-optimal paths evaporates with time while the pheromone on the near-optimal/optimal paths is reinforced, thus 
influencing more ants to follow the path [18]. 

In ACO, a population of agents (artificial ants) incrementally constructs solutions to combinatorial optimization 
problem by traversing a graph that encodes the optimization problem. Ant Colony system (ACS) algorithm, introduced by 
Dorigo and Gambardella [19] had been successfully deployed for discrete combinatorial optimization problems such as 
routing, and clustering. ACO algorithm has been extended to solving continuous combinatorial optimization problems 
using a variety of ACO algorithm called Variable Sampling Ant Colony Optimization (SamACO) algorithm [20]. 
SamACO algorithm offers an efficient incremental solution construction method based on the sampled values. 

The basic idea behind SamACO algorithm is that a population of agents (artificial ants) incrementally constructs 
solution to the sampled combinatorial optimization problem. The construction phase is guided by heuristic information and 
existing pheromone, which holds information about parts of a solution that have led to good results in the previous 
generation of ants. By means of pheromone update and transition rules, ACO probabilistically select the components 
values to concentrate the search in the regions of high quality solutions. 

The steps that are involved in the implementation of SamACO algorithm are as follows: 

Initialization Step 

The search space is bounded such that the decision variables (solution components) X ; has 
values Xj e [/, ,«,],/ = 1,2 , . . . , S , where /, and u j are the lower and upper bounds of the decision variables X, respectively, and 

S is the number of decision variables. The initial values of the decision variables are randomly sampled in the feasible 
domain as follows: 



in the neighborhood of the best-so-far solution per iteration, (m + 19) is the initial number of candidate values for each 
decision variable xj , rand is a random number uniformly distributed within [o, l] , i = 1, 2 S and j = 1, 2, . . . , (m + 1 ?) . 



foraging behavior of natural ants. It is a probabilistic technique for solving combinatorial optimization problems that can 
be reduced to finding good path through graphs. In their search for food, natural ants leave a trail of chemical substance 
called pheromone on the path they traverse in order to guide future ants toward optimal paths to food. The pheromone on 



x.^ = /, -1 — 1 -(j - 1 + rand. 4) 



m + $ 



(14) 



Where m is the number of ants, z? is the exploitation frequency which controls the number of values to be sampled 
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For each decision variable X , there are k. j sampled values xf\ x (2) , . .., x' k ‘ ’ from the continuous domain [/, , u , ] . 
Each solution component has associated a pheromone value, r/ , i = \,2,...,S , j = l,2,...,k i , and a component-pheromone 
matrix M can be generated. The pheromone value r/ reflects the desirability of adding the component value x‘ to the 
solution 



M = 



'{«} ' 

kvi fe 2i .4 8 }...k 8 ,k i ) 



(15) 



Transition 

The transition of the ants from one position to another is partially probabilistic and partially deterministic. An 
artificial ant k has a memory of the positions that it has already visited and the pheromone content at each location stored 

in a Tabu list T k . The memory size of the Tabu list depends on the ant population size as well as the number of movement 
made by the ants. In general, if there are m ants making N movement, the size of the Tabu list is ( mxN ). The iteration 
index Z* of the variable value selected by ant k for the i th variable is: 



1 if q < q o and 



l = 



• 0 if q^q o and 

f-f * if q>q 0 



j = f 



j * r 



(16) 



Where 

7 * = arg maxjr’ 1 *, t\ 2 \ ■■■ , (17) 

i = l,2,...,,s, k = 1, qe [o, l] is a uniform random value, and e [o, l] is a threshold parameter that 

represents the relative preference for either exploitation or exploration. 

Dynamic Exploitation 

When q<q„, the ant chooses exploitation in the neighborhood of the solution set with the highest pheromone 
value from the m solutions generated in the previous iteration. The dynamic exploitation is used as a local search method 
to fine-tune the best-so-far solution. A radius r confines the search in the neighborhood of the best-so-far solution 

x * 0) = (x[ o) , x°, x° s ) to the interval [x!°* -r n x® + r;.], i = l,2,--,S . The values of the variables in the best-so-far solution 
set are randomly selected to be increased, unchanged or reduced as 



r t .CT i ,u i ), 


0 < q < - 
3 






1 ^ 2 


(18) 




3 3 




■ r, .<7 ( , f), 


— < q < 1 
3 
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Where <7 e [ 0 , l] . 



The best-so-far solution set is then updated using elitism and generational replacement. The new solution set 
x = (xj, x 2 , x s ) is evaluated and used to replace the best-so-far solution set if there is an improvement in the result. The 
dynamic exploitation process is repeated for 7? times, and the newly generated solution components are recorded as 
x\'\ where j = m + l, m + 2, m+ g t , / = 1,2,. ..,5. The number of solution components generated during the dynamic 
exploitation process is denoted by g. . The radii are adaptively extended or reduced based on the exploitation result. If the 
best solution set produced by the exploitation process is better than the prevailing best-so-far solution set, the radii will be 
extended. Otherwise, the radii will be reduced. 



f r i .v e , v e >l 

j r. . v r , 0<v r <l 



09) 



Where v e and v r are the radius extension rate and the radius reduction rate, respectively. The initial radius value is 



given by: 



u. — /. 
2m 



( 20 ) 



Random Exploration 

When q > q o , the ant resort to probabilistic exploration to select a random index /)/' e {0,l,---,m+g ( }. Moving from 
position i , the ant chooses its next position j among the positions that have not been visited yet according to the 
probability distribution given as follows: 



Jj) _ . 






j = 0, 1, • ■ • , m + j 



( 21 ) 



The solution components of the worst X solution sets that are constructed by the ants in the previous iteration are 
discarded and each solution component is replaced by X new values generated by a random exploration process. If the 
worst solution sets are denoted by jc*'””" 1 *, x^ . The new solution components for the solution set x^are 

randomly generated as follows: 

x^ = /, + ( 77 ,. - 1 j). rand ■ j) (22) 



Where i = 1, 2, . . . , S and j = (m - x + 1), (m — x + 2), m . 



Random exploration ensures diversity and prevents premature convergence to local minima. 

Pheromone Update 

Initially, each solution component is assigned an initial pheromone value t o . The pheromone values are updated 
based on the quality of solutions constructed by the ants. The update is biased towards the best solutions constructed by the 
ants such that ACO concentrates the search in the regions of high quality solutions. Similar to MAX -MIN ant system [21], 
the pheromone values in SamACO algorithm are bounded to the interval [7 mm , 7’ ilia< ] ; in this case [ 0 . 1, l] . 
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The pheromones on the non-optimal paths are evaporated. The selected non-optimal solution components have 
their pheromones evaporated as 

7 , 0) <- (l - py/' + PT™ for 0<P<1 

i = 1, 2 , S and j = \,2,...,m (23) 

Where r mm is the predefined minimum pheromone value and p is the pheromone evaporation rate. 

The pheromones on the near optimal paths are reinforced, thus influencing more ants to follow the paths and 
hopefully find better solutions. The solution components in the selected best V P solutions have their pheromone reinforced 
as 

*■,")<- (1 - 0K 1 ' + for 0 < /? < 1 (24) 



i =1,2,. ..,5 and j = 1,2,..., 'T 

Where t ]m is the predefined maximum pheromone value, />’ is the pheromone reinforcement rate, and *P is the 
elitist number. 

In each iteration, the pheromone values of the solution components of the iteration-best are updated as the 
solution components of the best-so-far solution after the fitness evaluation of the m solutions constructed by the ants. 

IMPLEMENTATION 



Using MATLAB software, the proposed SamACO algorithm was implemented to compute the optimal switching 
angles that eliminate 5 th , 7 th , 1 1 th , and 13 th harmonics in an 1 1-level inverter. In this work, the population size is 40, and the 
number of iterations is 100. The solutions were computed by incrementing the modulation index, m j in steps of 0.001 from 
0 to 1. A personal computer (2.66 GHz Intel Core i7 processor with 4GB Random Access Memory) running MATLAB 
R2014b on OS X Yosemite version 10.10 was used to carry out the computations. 



The solution set at each step is evaluated with the fitness function. The objective here is to determine the 
switching angles such that the selected low order harmonics are either eliminated or minimized to an acceptable level while 
the fundamental voltage is obtained at a desired value. For each solution set, the fitness function is calculated as follows 
[ 12 ]: 



/ = ] 



17 v:-vP 


s i ( v, V 


100 1 


+ 1 — 50 -^- 


[1 u 


& h , { v, ) \ 



1 = 1,2,. ...S 



(25) 



Subject to 

0 < a, < a 2 < ... < a s < 



Where V," is the desired fundamental output voltage, S is the number of switching angles, h s is the order of the s"‘ 
viable harmonic at the output of a three phase multilevel converter. For example, h 2 = 5 , h t = 11. It should be noted that 
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different weight are assigned to different harmonics in eqn. (25). Each harmonic ratio is weighted by inverse of its 
harmonic order, i.e. l/h s . By this weighting method, higher importance is assigned to the low order harmonics, which are 
more harmful and difficulty to remove with filter. 

In order to validate the observed analytical results, an 11 -level single-phase Cascaded H-Bridge inverter was 
modelled in MATLAB-SIMULINK using SimPower System block set. In each of the five H-Bridges in the 11-level 
single-phase Cascaded H-Bridge inverter, 12V dc source is the SDCS, and the switching device used is Insulated Gate 
Bipolar Transistor (IGBT). Simulations were performed at the fundamental frequency of 50 Hz using the solution set found 
at the modulation index, m t of 0.795: cr, = 6.99" , a, = 1 9.05“ , a, = 28.0 1" , a, = 45.99" , and a., = 62.61" . Fast Fourier 
Transform (FFT) analysis of the simulated phase voltage waveforms was done using the FFT block to show the harmonic 
spectrum of the synthesized AC voltage. 

RESULTS 

The plots of fitness value for each set of switching angles versus modulation indices over the range of 0.1 to 1.0 is 
shown in Figure 3. For values of the fitness function less or equal to 0.01, SHE equations are solvable, otherwise they are 
unsol vable. 




Modulation Index 

Figure 3: Fitness Function at Various Modulation Indices 

As shown in Fig. 4, there are multiple solution sets at some modulation indices. In such cases, the solution set 
with the least THD value is chosen. 




0.4 0.6 

Modulation Index 



+ alphal 
alpha2 - 
+ alpha3 
alpha4 
alpha5 ' 



Figure 4: Switching Angles at Various Modulation Indices 

At low modulation indices [0, 0.256], high modulation indices [0.934, 1] as well as some modulation indices, no 
solution sets are available. For those modulation indices, it is either there is no solution set or SamACO could not find one. 
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The former reason is more plausible than the latter. 




Figure 5: THD Values at Various Modulation Indices 

From the simulated output voltage waveform, the peak value of the fundamental output voltage is 60.74T (/) u() . This 
value is in close agreement with the analytically computed value given by eqn. (8) as; 



V. = m. 



4 sV* 



K 



= 0.795 



4x5x12 



K 



= 60.73V, 



( peak ) 



The FFT analysis of the synthesized voltage waveform shown in Fig. 6 reveals the complete elimination of the 
5"\ 7'*, \Y h and 13" 1 harmonics as their values tend towards zero. 



The THD in line-to-line voltage as computed analytically and from simulation are 4.54% and 4.58% respectively. 
The analytical and simulation values of THD are in close agreement thereby validating the analytical results. The THD 
value of 8.01% shown in Fig. 6 is the THD value of the synthesized voltage waveform of a single phase inverter which 
includes triplen harmonic components. 




Figure 6: (a) Synthesized Waveform (b) Harmonic Spectrum of an 11-Level Inverter with m i = 0.795 



CONCLUSIONS 

Variable Sampling Ant Colony Optimization (SamACO) algorithm with random initial values has been 
successfully implemented for solving the transcendental nonlinear equations characterizing the harmonics in an 11-level 
inverter. The proposed method is derivative-free, accurate and globally convergent. The absence of the selected harmonics 
in the synthesized output phase voltage validates the analytically computed results. 
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